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^ ' Abstract 

' Azimuthal angle decorrelation in inclusive dijet cross sections is stud- 

, ied analytically to take into account the next-to-leading corrections to the 

' BFKL kernel while keeping the jet vertices at leading order. The spectral 

representation on the basis of leading order eigenfunctions is generalized to 
\£) ■ include the dependence on conformal spins. With this procedure running 

coupling effects and angular dependences are both included. It is shown 
how the angular decorrelation for jets with a wide relative separation in 
■ rapidity largely decreases at this higher order in the resummation. 

Q ■ 1 Introduction 



Among the many relevant questions still open in Quantum Chromodynamics a 
very interesting one is how to describe scattering amplitudes in the so-called 
Regge limit. In Regge asymptotic the center-of-mass energy, s, is much larger 
than all other Mandelstam invariants and mass scales present in the process 
under investigation. It is possible to perturbatively keep track of the differ- 
ent contributions to the amplitude if some hard scale is present so that the 
strong coupling remains small. It is then needed to rcsum logarithmically en- 
hanced contributions of the form (a s Ins)™ to all orders. This is achieved using 
the leading-logarithmic (LL) Balistky-Fadin-Kuraev-Lipatov (BFKL) evolu- 
tion equation pQ . 

The BFKL approach predicts a power-like rise in s of total cross sections. A 
lot of attention has been given to the search for BFKL effects in deep inelastic 
scattering (DIS) due to the rapid growth of structure functions at small values of 
Bjorkcn x. However, the golden process where the resummation of Ins is most 
important is the total cross section of two photons with large and similar virtu- 
alities. In this configuration the small x resummation, which includes ordering 
in rapidity and not in transverse scales, should be the dominant contribution to 
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the scattering. This is not necessarily the case in DIS where ordering in kt is 
important given that parton evolution takes place between an object with large 
transverse size, the proton, and a small highly virtual photon. 

Observables where BFKL effects should prevail then require of enough energy 
to build up the parton evolution, and the presence of two large and similar 
transverse scales. In this work an example of this kind is investigated in detail 
by analytic means: the inclusive hadroproduction of two jets with large and 
similar transverse momenta and a big relative separation in rapidity, the so- 
called Mueller-Navelet jets. When Y, the distance in rapidity between the 
most forward and backward jets, is not large a fixed order perturbative analysis 
should be enough to describe the cross section but when it increases a BFKL 
resummation of (a s Y) n terms is needed. 

Mueller-Navelet jets were first proposed in Ref. [5] as a clean configuration 
to look for BFKL effects at hadron colliders. A typical power-like rise for 
the partonic cross section was predicted in agreement with the value of the 
asymptotic LL hard Pomeron intercept. However, at hadronic level, forward and 
backward jets are produced in a region of fast falling of the parton distributions, 
reducing the rise of the cross section. A way to make small x resummation 
effects more explicit is to look into the azimuthal angle decorrelation of the pair 
of jets. The relevant subprocess is parton + parton — ► jet + jet + any number 
of soft emissions inside the rapidity interval separating the two jets. BFKL 
enhances soft real emission as Y increases reducing in this way the amount of 
angular correlation originally present in the back-to-back in transverse plane 
Born configuration. The LL prediction for this azimuthal dependence was first 
investigated in Ref. 0. 

The results at LL are known to overestimate the rate of decorrelation and to 
lie quite far from the experimental data 0] as obtained from the Tevatron and 
subleading higher order effects have been called for an explanation of this dis- 
crepancy. Running coupling effects and kinematic constraints have been consid- 
ered in Ref. |Sj ■ In the present work the main target is to analytically understand 
how to include the a s (a s Y) n next-to-leading logarithmic (NLL) corrections to 
the BFKL kernel The effects of this kernel were numerically investigated 
using an implementation [7] of the NLL iterative solution proposed in Ref. 8 
(different reviews can be found in 9 ) . It is left for future analysis the inclusion 
of the next-to-leading order (NLO) jet vertex [Hi], the investigation of more 
convergent versions of the kernel ^5 an d a study of parton distributions effects 
in order to make reliable phenomenological predictions at a hadron collider. 
Mueller-Navelet jets should be an important test of our understanding of small 
x resummation to be performed at the Large Hadron Collider at CERN. 

After this brief Introduction, in Section 2 the normalization for the gluon 
Green's function is indicated together with the formulae for the partonic cross 
section. Then the operator formalism suggested by Ivanov and Papa in Ref. J2| 
is extended to introduce angular dependences. The form of the NLL kernel for all 
conformal spins calculated by Kotikov and Lipatov in Ref. ^2] is also discussed. 
Towards the end of the section a compact expression for the angular differential 
cross section which includes the NLL contributions is derived. In Section 3 the 
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numerical study of the previous formulae is discussed in detail. Finally, several 
Conclusions are drawn and different lines for future research highlighted. 



2 Calculation of the dijet partonic cross section 

As indicated in the Introduction, in this analysis the object of interest is the 
partonic cross section parton + parton — > jet + jet + soft emission, with the 
two jets having transverse momenta q\ and q 2 and being produced at a large 
relative rapidity separation Y. This can be related to the external hadrons by 
its approximate relation to the longitudinal momentum fractions carried by 
the jets, i.e., Y ~ In x \x 2 sj ygf gf, and a convolution with parton distribution 
functions whose analysis is left for a future work. In the present framework the 
resummed differential partonic cross section for the particular case of gluon- 
gluon scattering is 



da _ n 2 a 2 s f(q\,q 2 ,Y) 



d 2 qxd 2 q 2 2 q\q\ 



(1) 



with the usual notation a s — a s N c /iT. One can now introduce the Mellin 
transform of the BFKL gluon Green's function in rapidity space: 

/ (qi , Y) = J (<li , &) ■ (2) 

The normalization for the BFKL integral equation, including for simplicity only 
the LL terms, then reads 



ufu (qi, 92) = <5 (2) (<fi - q 2 ) 



dk I /-. \ qjU(qi,q2) I ( ^ 
I ^2 \U[k,q 2 )- z -2 . (3) 

7r(gl-fcJ ^ k 2 +[q l -k) ) 

As it is well known, including the angular dependence on the transverse plane 
of qi and q 2: the LL solution to Eq. J3Jl can be written as 

1 00 r°° -i-i P m(e l -e 2 ) 

= AE *»{<8) 2 (ql) w ~"—- — ju—, (4) 



n— — 00 



where the eigenvalue of the LL kernel, 

Xa{n,u) = 2^(l)-vQ + ^+^) , (5) 

is expressed in terms of the logarithmic derivative of the Euler Gamma function. 
As it stands the n variable corresponds to a Fourier transform in the angular 
sector. A more sophisticated interpretation arises if Eq. J2J is considered for 
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non-zero momentum transfer (see, e.g., |14|L In this case, if a representation in 
the complex plane for the transverse momenta of the form q = q x + iq y is intro- 
duced, it can be shown how the BFKL equation corresponds to a Schrodinger- 
like equation with a holomorphically separable Hamiltonian where — iY is the 
time variable. Both the holomorphic and antiholomorphic sectors in the Hamil- 
tonian are invariant under spin zero Mobius transformations with eigenfunctions 
carrying a conformal weight of the form 7 = ^ + iv + ^ . In the principal series 
of the unitary representation v is real and \n\ the integer conformal spin [15j . 

The partonic cross section is obtained by integration over the phase space 
of the two emitted gluons together with some general jet vertices, i.e. 

a{a s ,Y,p\ 2 ) = J d 2 qi J d 2 q 2 $ jctl (qi,P?) $jct 2 (g2,pl) d 2^p^ 2 • ( 6 ) 

In the perturbative expansion of these jet vertices, $j e t, = 5>jg t . + c*s^jet, + • • •) 
only leading-order terms are kept: 

= 0{q 2 -p 2 ), (7) 

where p 2 corresponds to a resolution scale for the transverse momentum of the 
gluon jet. In this way a full NLO accuracy is not achieved but it is possible to 
pin down those effects stemming from the gluon Green's function. To extend 
this analysis it would be needed to calculate the Mellin transform of the NLO 
jet vertices in Ref. |1U) where the definition of a jet is much more involved than 
here. Therefore one can proceed and write 

( g£ (gi.Pj) gjj&Plh - . y"1 (XS 

a(a s ,Y,p 12 ) = —— d qx d q 2 — 5 2 / («i> «2> Y ) • (8) 

^ J J <7i 92 

At this stage it is very convenient to recall the work of Ref. ^21 an d to introduce 
the following transverse momenta operator representation: 

Q Wi) = Qi \Qi) (9) 

with the normalization 

(qi\lW 2 ) = 5^(qi-q 2 ). (10) 
In this notation the BFKL equation simply reads 

(u-K)f u = 1 (11) 



where the kernel has the expansion 

K = a s K + a 2 s k x + ... (12) 
To NLO accuracy this implies that the solution can be written as 

L = (w - u s kv) 1 + a 2 s (u - a s k ) 1 Ki (u - a s Ko) ' + O (a 3 s ) . (13) 
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The next step is to define a basis on which to express the cross section. To 
generalize the study in Ref. J^l this basis should carry not only the dependence 
on the modulus of the transverse momenta but also the dependence on their 
angle on the transverse plane: 



(q\ v, n) 



1 



irV2 



2\ w ~2 e in6 



(14) 



The projection (n, v \q) would be the complex conjugate of the previous expres- 
sion. This basis has been chosen such that it is orthonormal: 



(n', v \v, n) 



5 {y — v') 5„ 



(15) 



The action of the NLO kernel on this basis, which was calculated in Ref. |13| . 
contains non-diagonal terms and can be written as 



K\u,n) = (a s Xo(H,v)+a 2 s xi{\n\,p) 



■a. 



d 



2xo (N ,u) i— +log^ 2 + i- 5 -Xo(\n\,v) 



t) 



d 



W,n), (16) 



where, from now on, a s stands for a s (/i 2 ) , the coupling evaluated at the renor- 
malization point \i in the MS scheme. The first line of Eq. I|16[l corresponds to 
the scale invariant sector of the kernel. The function xi for a general conformal 
spin reads 

Xi(n,7) = 5xo(^,7) + ^C(3)-^-Xo(«:7) 



r (7 + 

7r 2 cos (wy) 



I-7 + -) -20(n, 7 )-20(ra,l- 7 ) 



4 sin 2 (tt7) (1 - 2j) 



- 1 



n,\ 2 + 37(1-7) 
" NlJ (3 -2 7 ) (l + 2 7 ) 

7(1-7) £ 



>n() 



N?J 2 (3 - 2 7 ) (1 + 2 7 ) 



(17) 



The definitions S = (4 - tt 2 + 5/3 a /N c ) /12, and /3 = (11N C - 2n/)/3, have 
been used. The function <f> can be found in Ref. |13| and reads 



(-l)(fc+!) 



ip'(k + n+l)-i/j'(k + l) 



+ n + 1) + + 1)) + ^ + l)-^ + n + l) >| 



(18) 



with 



4/3' (7) 



1 + 7 



(19) 
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In this basis terms with derivatives are associated to the running of the 
coupling ^B], this is the case in the second line of Eq. i|16|) . In particular, the 
part containing 2^7, which breaks the v — > — v symmetry, will be shown to give 
a zero contribution to the cross section when p\ = p\. It will also be shown how 
the term with «Xo^7 mixes in a non-trivial way the Green's function with the 
jet vertices. 

To represent the cross section in the present formalism the starting point is 
to project the jet vertices on the basis in Eq. (fTl)l : 

d2 (&pj) m n) = i i_ s ci (20) 

<r V2 (| - iv) 

The C2 (v) projection of ^jet 2 on ( n i v \ Q) is the complex conjugate of l|2U|) with 
p\ being replaced by p\. The corresponding inverse relations are 

J6tl 2 = £ / ^CiMtoKH?), (21) 



q 2 

n— — oo 



£ / ^c 2 (^)5 n , (g \v,n). (22) 

— J OO 



The cross section can then be rewritten as 



^2^2 00 roo f-oo 

a (a s ,Y,p\ 2 ) = — — ^ ^2 diy dv'ci(i>)c2(i' / )8 ni o5 n >$ 

* . /_ J — oo J —oo 

■e" Y (n,H fl\u',n'}. (23) 



n,n — — oo 



2ni 



Making use of the operator representation in Eq. I|13|) , the action of the kernel 
in Eq. I|l(j|) and integration by parts, a can be expressed as 

2-2 °° poo 

a(a s ,Y,p 2 12 ) = — ^ / ^e«»(H.)Y CiMc2(l/) ^ (24) 



n— — oo 



Xl(l „ M + ^( log(M2) + i| log (|M) + «|) M(l „ u) 

For the LO jet vertices the logarithmic derivative in Eq. (|24|l explicitly reads 



x i 1 + a? Y 



The angular differential cross section can be calculated considering the following 
representation of the C1C2 product: 

*<■"«<■"*- - ^mwUiS)" (26) 



G 



with tf> = Q\ — 62 — I"- Therefore, in the case where the two resolution momenta 
are equal, p\ = p 2 = P 2 , the angular differential cross section can be expressed 
as 



da (a S: Y,p 2 ) _ ir 2 g 2 s ^ J_ . 
cLf> ~ 4p 2 2tt 6 



n— — 00 



dv e &.Xo{\n\, V )Y 



l + a 2 Y 



A X0 (| n |, I/ )[ 2 log(^ 




The term proportional to in Eq. I)24|l gives no contribution after integra- 
tion as it is an odd function in v. Within NLO accuracy there is freedom to 
exponentiate the integrand of this result. In this work this is done both for the 
scale invariant and for the scale dependent terms, in close resemblance with the 
property of reggeization. Furthermore, renormalization group improved pertur- 
bation theory is called for to introduce the replacement 



a -wM^ ^ & - {p2) - (28) 



In this way the differential distribution can be conveniently rewritten as 

e'^C n (Y) , (29) 



da (a s ,Y,p 2 ) _ ^g 2 1 

cLfr ~ 2p 2 2^ ^ 



n— — 00 



with 



a 3 (p 2 )Y^ X o(|n|^)+a s (p 2 )^Xi(|r l |^)-^i!p^ 



Cn (Y) = r d £ e - : rr — — . (30) 

Different interesting observables can be constructed with these coefficients and 
they will be studied in detail in the next section. Due to the large and negative, 
with respect to the LL terms, size of the NLL corrections it will turn out that 
the exponentiated form in Eq. (|30l) is mandatory in order to reach convergent 
results. This is discussed below. 



3 Analysis of convergence and study of observ- 
ables 

To investigate the properties of the different parts of the kernel it is useful to 
introduce five types of coefficients. The first one is 

f°° dv p««Yxo(|n|,iO 

cm - /JW. m 
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defined at LL accuracy and previously studied in the literature (see last refer- 
ence in [3|). When all the NLL terms are exponentiated as in Eq. (|30|l it will 
be referred to as C^ LL , while if the NLL pieces are not exponentiated, as in 
Eq. (SJ, it is named C^ xpansion . If in Eq. (JSJ the xi kernel is removed then 
the coefficients correspond to the case of LL plus running coupling and it is 
noted as cR unmn g cou P lm s. Finally, the scale invariant contributions can be iso- 
lated by setting (3q to zero in the exponent of Eq. 1|30[) and this will be called 

/''Scale invariant 

The coefficient governing the energy dependence of the cross section corre- 
sponds to n = 0: 

a(a s ,Y,p 2 ) = ^rCo(Y). (32) 

As the NLL corrections are large and negative, when not exponentiated they 
lead to a non-convergent behavior. This is seen in Fig. ^ where the resolution 
scale p — 30 GeV has been chosen. The values rif — 4 and Aqcd = 0.1416 GeV 
were taken in a s (p 2 ) = 4N C / (/3 In (p 2 /Aq CD )) ■ In this plot it can be seen how 
the coefficient £^ x P anslon generates an unphysical behaviour in contrast to the 
exponential rise associated to the other coefficients. Nevertheless it should be 
noticed that the term proportional to a 2 s Y 



(33) 



would be added to xi in Eq- <P7|) if the NLO jet vertex +a s c^ was brought 
into the calculation. The convergence properties of Eq. i|27[) might well change 
in that case. 

A familiar consequence of introducing the effects of the running of the cou- 
pling is that the LL intercept is reduced as can be seen also in Fig.^ Meanwhile, 
if the scale invariant sector of the NLL kernel, Xii is a ls° introduced then a fur- 
ther decrease of this rise takes place. The n = coefficient is directly related to 
the normalized cross section 

£00 CoOO 

*(0) Co(0)' [ ' 




The rise in rapidity of this observable is shown in Fig. [21 Clearly the NLL 
intercept is very much reduced with respect to the LL case. In Fig. |31 large 
values of Y are considered in order to approach the asymptotic regime. The LL 
intercept tends to the asymptotic value of 4 a s (30) In 2 ~ 0.37 while the NLL 
result lies around one third of this number. 

Independently of the NLL kernel being exponentiated or not, the remaining 
coefficients with n > 1 all decrease with Y. This can be seen in the plots of 
Fig. The consequence of this decrease is that the angular correlations also 
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Figure 3: Intercepts for the cross section as a function of rapidity. 



diminish as the rapidity interval between the jets gets larger. This point can be 
studied in detail using the mean values 



(cos (mcf))) = 



C m (Y) 
Co (Y) ' 



(35) 



(cos ((f))) is calculated in Fig. EI The most important consequence of this plot is 
that the NLL effects dramatically decrease the azimuthal angle decorrelation. 
This is already the case when only the running of the coupling is introduced but 
the scale invariant terms make this effect much bigger. This is encouraging from 
the phenomenological point of view given that the data at the Tevatron typi- 
cally have lower decorrelation than predicted by LL BFKL or LL with running 
coupling. It is worth noting that the difference in the prediction for decorrela- 
tion between LL and NLL is mostly driven by the n = conformal spin. This 
can be understood looking at the ratio 



NLL 



(cos (</>)) 
(cos((/>)) LL 



C NLL (y) C LL ( Y ) 

C£ LL (Y) C : LL (Y) ' 



and noticing that 



ct^y) 

Li > (Y) > L 



(36) 



(37) 



This ratio is calculated in Fig. El This point is a consequence of the good 
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Figure 5: Dijet azimuthal angle decorrelation as a function of their separation 
in rapidity. 
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Figure 6: Comparative ratio between the NLL and LL coefficients for n = 1 
conformal spin. 
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convergence in terms of asymptotic intercepts of the NLL BFKL calculation for 
conformal spins larger than zero. In particular the n — 1 case is special in that 
the property of zero intercept at LL, xo(L 1/2) = 0, is preserved under radiative 
corrections since 

is also zero. For completeness the m = 2,3 cases for (cos(m^)) are shown in 
Fig. [7| These distributions are relevant because they prove the structure of 
the higher conformal spins. The trend is the same as previously discussed: the 
correlation increases as higher order corrections in the small x resummation are 
included. 

4 Conclusions 

An analytic procedure has been presented to calculate the effect of higher order 
corrections in the description of Mueller-Navelet jets where two jets with mod- 
erately high and similar transverse momentum are produced at a large relative 
rapidity separation in hadron-hadron collisions. This is a promising observ- 
able to study small x physics at the Large Hadron Collider at CERN given its 
large energy range. The focus of the analysis has been on those effects with 
direct origin in the NLO BFKL kernel, while the jet vertices have been con- 
sidered at LO accuracy. It has been shown how the growth with energy of the 
cross section is reduced when going from a LL to a NLL approximation, and 
how the azimuthal angle decorrelations largely decrease due to the higher order 
effects. The present study has been performed at partonic level while the im- 
plementation of a full analysis, including parton distribution functions, NLO jet 
vertices and the investigation of collinearly improved kernels, will be published 
elsewhere. 
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